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ABSTRAGT 

We have used the publicly released Dark Energy Survey data to hunt for new satellites of the Milky 
Way in the Southern hemisphere. Our search yielded a large number of promising candidates. In this 
paper, we announce the discovery of 9 new unambiguous ultra-faint objects, whose authenticity can be 
established with the DES data alone. Based on the morphological properties, three of the new satellites 
are dwarf galaxies, one of which is located at the very outskirts of the Milky Way, at a distance of 
380 kpc. The remaining 6 objects have sizes and luminosities comparable to the Segue I satellite and 
can not be classified straightforwardly without follow-up spectroscopic observations. The satellites 
we have discovered cluster around the LMG and the SMG. We show that such spatial distribution 
is unlikely under the assumption of isotropy, and, therefore, conclude that at least some of the new 
satellites must have been associated with the Magellanic Glouds in the past. 

Keywords: Galaxy: halo, galaxies: dwarf, globular clusters: general, galaxies: kinematics and dynam¬ 
ics 


I. INTRODUCTION 

Gurrently, there are no strong theoretical predictions 
as to the luminosity function and the spatial distribution 
of the Galactic dwarf companions. In other words, to¬ 
day’s semi-analytic models are so flexible they can easily 
produce any number of Milky Way (MW) satellites: from 
a few tens up to several thousands (e.g. iKoposov et al.l 
|2QQ9[ ). Similarly, a large range of spatial arrangements 
is possible: from nearly is otropic to strongly planar (e.g. 
Ib^ V Baumgardtl[2Qd^ . This is why to improve our 
understanding of the physics of the Universal structure 
formation on small-scales we look to observations. How¬ 
ever, to date, only a third of the sky has been inspected. 
The Sloan Digital Sky Survey (SDSS) has demonstrated 
the power of deep wide-area imaging to fuel resolved stel¬ 
lar populations studies of the Galactic halo. The analysis 
of the SDSS data covering only 1/5 of celestial sphere has 
more than doubled the number of the known Galacti c 
dwarf satellites (see e.g. IWillmanll2QIQl: lBelokur^l2QI3[ ). 
Recently, both PanSTARRS and VST ATLAS have ex¬ 
panded the surveyed region signific antly, but found little 
of the wealth they c ame to seek (jLaevens et al.l l2Ql4 
iBelokurov et al.l[2QI4f ). 

The SDSS discoveries extended the dwarf galaxy 
regime to extremely low luminosities and sizes. As a 
result of adding these new faint surface-brightness ob¬ 
jects (known as ultra-faint dwarfs, UFDs) to the panoply 
of Milky Way companions, it has been revealed that 
there appears to be a gap in the distribution of effective 
radii between globular clusters and dwarfs which extends 
across a large range of luminosities. In the absence of a 
working definition of a galaxy, an ad hoc combination of 
morphological, kinematic and chemical properties is re¬ 
quired to be classified as one. Therefore, some objects 
that live dangerously close to the notional gap may risk 
misclassification. These faintest of the UFDs are only 
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detectable with surveys like SDSS out to a small frac¬ 
tion of the MW virial radius. As a result, their total 
number in the Galaxy must be reconstructed under the 
assumption of the shape of their galactocentric radial dis¬ 
tribution. Given that only a handful of such objects are 
known today, their radial profile is not observationally 
constrained. However, if it is assumed that the faintest 
of the UFDs are linked to the small-mass field dark mat¬ 
ter (DM) sub-halos, their distribution can be gleaned 
from cosmological zoom simulations. This is an example 
of so-called sub-halo abundance matching which predicts 
that the bulk of the Galactic dwarf galaxy population is 
in objects with luminosities My < —5. The details of 
the semi-analytic calculations may vary but the result 
remains most astounding: there ought to be hundreds 
of galaxies with My ^ —I, if the faintest of the UFDs 
occupy field DM su b-halos (see e.g. iTollerud et al.l[200^ 
iBullock et al.ll2QTQ[ ). 

It is, however, possible that the smallest of the UFDs 
may have been acquired via a different route. It has been 
suggested that some of these objec ts may have been ac¬ 
creted as part of a group (see e.g. iBelokuro^l2QI3f ) . In 
the case of the UFD Segue 2 for example, the group’s 
central object appears to have been destroyed and can 
now be detected only as tidal debris in t he MW halo 
(jBelokurov et al.l[2QQ9l: iDeason et al.ll2Ql3) . In this sce¬ 
nario, the satellites of satellites survive (unlike their pre¬ 
vious host) in the MW halo, but their spatial distribu¬ 
tion should differ dramatically from that predicted for 
the accreted field sub-halos. Not only the radial density 
profile is modified, the satellites’ morphological proper¬ 
ties might have been sc ulpted by the pre-processing (see 
e.g. IWetzel et al.l[2QI5[) . Thus, the total numbers of the 
UFDs would be biased high if estimated through sim¬ 
ple abundance matching calculation, while their spatial 
anisotropy would be greatly under-estimated. 

On the sky, the remaining uncharted territory lies be¬ 
neath declination S = —30° and is currently being ex- 
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Figure 1. SExtractor star/galaxy indicator SPREAD .MODEL 
as a function of the PSF magnitude. The red dashed curves 
are the lines of constant |SPREAD_MODEL| =0.003, while the 
red dotted lines are delineating the region used to select stars 
|SPREAD_MODEL| <0.003+SPREADERR_MODEL 

plored by the rapidly progressing Dark Energy Survey 
(DES). This is the realm of the Magellanic Clouds, long 
suspected to play a role in s haping the Galacti c satel¬ 
lite population. For example, iLvnden-Belll (|1976[ ) points 
out the curious alignment between several of the Milky 
Way dwarfs and the gaseous strea n i ema nating from the 
Magellanic Clouds. iKroupa et al.l (|2QQ5[) developed this 
picture, adding more Milky Way satellites to the Great 
Magellanic family. The idea also receives support from 
the planar structures of sa tellites around M31 and other 
nearb y galaxies (see e.g., Ilbata et ^ 120131: iTnllv et al.l 
|2Q15[) . To explain the peculiar orbital arrangement of the 
Large Magellanic Cloud (LM C) and the fainter dwarfs 
within the ACDM paradigm, iD’Onghia fc Lake! (|2QQ8f ) 
consider a cosmological zoom simulation of accretion 
of an LMC-type object onto a Milky Way-size galaxy. 
In their simulation, approximately 1/3 of the satellites 
around the prototype MW at redshift z = 0 come from 
the LMC group. They conjecture that the process of 
dwarf galaxy formation at the epochs around the re¬ 
ionisation is favoured in group environments, suggesting 
that at least in part, this can be explained by a signifi¬ 
cant difference in gas cooling times in the halo of a large 
galaxy like MW and a small group of the size of the LMC. 

In this paper, we announce the discovery of a large 
number of new faint satellites detected via an automated 
stellar over-density search in the DES public release data. 
The paper is structured as follows. Section 2 outlines 
the steps required to produce object catalogues from the 
calibrated DES image frames. Section 3 presents the 
details of the satellite search routine. Section 4 gives an 
overview of the properties of the detected satellites, while 
Section 5 discusses their spatial distribution. The paper 
closes with Discussion and Conclusions in Section 6. 

2. DES DATA PROCESSING 



Figure 2. Quality of the star-galaxy separation as inferred using 
the SPREAD_MODEL parameter. The photometry extracted from 
the DES images is compared to the star-galaxy separation in the 
CFHT Wide survey of the same area on the sky. Black solid curves 
show the fraction of objects classified as likely stars by CFHT 
CLASS-STAR >0.5 that are also classified as stars according to 
the SPREAD-MODEL cut (Eq. [T}. This is a measure of complete¬ 
ness of our stellar sample. Red dashed curve shows the fraction of 
objects classified as stars according to the SPREAD-MODEL cut 
but are classified as likely galaxies by CFHT CLASS-STAR < 0.5. 
This is a measure of the contamination. 

DES enjoys the use of the 2.2-degree field-of-view of 
the 570 megapixel DECam camera mounted on the 4m 
Blanco telescope at Cerro Tololo in Chile to obtain deep 
images of an extensive swath of the Southern sky around 
the Magellanic Clouds. DES is an international collab¬ 
oration, whose primary focus is on the extra-galactic 
science. DES started taking data in August 2013 and 
will continue for five years. Currently, all its data prod¬ 
ucts are proprietary, but the individual processed images 
are eventually released to the community through the 
NO AO infrastructure with a 1 year delay. In this paper, 
we present the results of the analysis of the entire DES 
imaging data released to date, i.e. ^2100 square degrees 
obtained in the year of operations. Here, we concen¬ 
trate on what clearly is a by-product of the DES survey: 
the stellar photometry and astrometry. 

The NO AO Science Archival! hosts both raw and pro¬ 
cessed DES images. We have queried the database 
through the interface provided and selected all avail¬ 
able InstCal calibrated images, the corresponding weight 
maps and the maskfl The images were first staged for 
the download, and subsequently copied across for fur¬ 
ther processing. The total size of the released images, 

^ http://WWW.portal-nvo.noao.edu/ 

^ The exact query used to select images from the NOAO Sci¬ 
ence archive: SELECT reference, dtpropid, surveyid, release-date, 
start-date, datc-obs, dtpi, ra, dec, telescope, instrument, filter, ex¬ 
posure, obstype, obsmode, proctype, prodtype, seeing, depth, dtac- 
qnam, reference AS archive-file, filesize, md5sum FROM voi.siap 
WHERE release-date < ’2015-01-18’ AND (proctype = ’InstCal’) 
AND (prodtype IS NULL OR prodtype <> ’png’)) AND (dtpropid 
ILIKE ’%2012B-0001%’ OR surveyid ILIKE ’%2012B-0001%’) 
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Figure 3. False-colour map of the density of the main sequence turn-off (MSTO) stars with 0.2<g—r<0.6 as observed by the DES survey. 
The density in different magnitude bins corresponding to different distances is used to create different colour channels: blue 17<r<19 
(corresponding to distances of ~ 4—lOkpc), green 19<r<21 (10—25kpc), red 21<r<22.75 (25—56kpc). The MW satellites discovered in 
this paper are shown with white unfilled circles. Note that the declination axis has been broken to avoid showing large empty area between 
-35< S <-5 without DES data 


weight maps and masks in Rice compressed format was 
^4.6 TB. Once uncompressed, the size of the images and 
weight maps increases by a factor of ^ 10. For the anal¬ 
ysis reported below, only g^r and i band images were 
analysed 30 TB of data). 

DES images are taken with 90s exposure and the DE- 
Cam pixel scale is 0.26"/pixel. InstCal products, de¬ 
livered by the DECam pipeline, are single-frame images 
(not stacks) that have been bias, dark and flat-field cal¬ 
ibrated as well as cross-talk corrected. Additionally, de¬ 
fects and cosmic ra^ have been masked and the WCS 
astrometry providedlj. 

2.1. From Images to Catalogues 

Eor the subsequent catalogue creation, we have 
relied heav ily on the SExtract or /PSEEx rou¬ 
tines (see iBertin fc Arnoritsl 119961: iBertinI 120111: 
lAnnnnziatella et al. I 120131 for more details). Typi¬ 
cally, processing a single frame consisting of 60 CCD 
chips with SExtractor/PSEEx took ^ 2.5 CPU-hours, 
leading to a total budget for the three DES filters of 
^15000 CPU-hours. The image processing has been 
carried out using the Darwin HPC computing facilit}^ 
as well as a local 8 node/96 core computer cluster. The 
main catalogue assembly steps are as follows: 

• Initial SExtractor pass. To provide the starting 
parameters for the PSEEx routine, the DES im¬ 
ages are analysed by SExtr actor. The flag settings 
are standard, albeit the PH0T_APERTURES flag is ad¬ 
justed to match the value of the G_SEEING keyword 
in the header (if it is available). Additionally, the 
downloaded weight maps are provided to SExtrac- 
tor via the WEIGHT_MAP option. 

• PSEEx run. The output of SExtractor is used to 
calculate the PSEs for the subsequent photometry. 
PSEEx is run on each CCD chip separately with 
default PSEEx options, except the PSF_SIZE key¬ 
word which is set to (40,40) pixels. 

^ Detailed information available through the DECam Data 
Handbook at http://bit.ly/lD7ZGDc 

^ http://www.hpc.cam.ac.uk 


• Einal SExtractor pass. Using the PSE models com¬ 
puted in the previous step, SExtractor is run again 
to determine Model and PSF magnitudes. This 
step is the longest in the entire sequence. 

• Ingestion. Eor each filter, the resulting SExtrac¬ 
tor catalogues are ingested into a PostgreSQL 
database as separate tables . Q3C spatial indices 
([Kodosov fc BartunovI 1200^ are created to speed 
up further steps. 

• Duplicate removal. Due to significant image over¬ 
laps, sources appearing more than once have to 
be flagged in the catalogues. Eor each source we 
perform a search with 1" radius, and mark all ob¬ 
jects coming from different frames and HDUs as 
secondary. 

• Photometric calibration. To convert instrumen¬ 
tal magnitudes into calibrated ones for each frame 
(including multiple HDUs) a cross-match between 
the DES and the APASS DR7 survey data is per¬ 
formed. The zero-point is then measured as a me¬ 
dian offset with respect to the APASS photome¬ 
try on a per-field basis. The resulting photometric 
precision as estimated using overlaps of different 
frames is: 0.03 mag (Gaussian a) in the g-band, 
0.025 mag in the r-band, and 0.033 mag in the i- 
band. 

• Band-merging. The lists of primary sources in the 

r and i bands are combined into the final cat¬ 
alogue using the matching radius of 1 arcseconds. 
Eor the purpose of this paper, the catalogue was 
based on r-band, i.e. r-band detection was required 
in order for the object to be in the final catalogue. 

• The final Q3C index is created on the band-merged 
catalogue to speed up spatial searches in the 
dataset. 

2.2. Photometry and star-galaxy separation 
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Table 1 

Parameters of the discovered MW satellites. 


Name 

a 


SigniP 

m-M’" 

Dist©’^ 

Mv 

Fmaj 

rl/2'' 

ri/2® 

ellipticity 

PA 

BE® 


[deg] 

[deg] 

icance 

[mag] 

[kpc] 

[mag] 

[arcmin] 

[arcmin] 

[pc] 


[deg] 


Reticulum 2 

53.9256 

-54.0492 

48.5 

17.4 

30 

-2.7T0.1 

o 07+0.23 
' -0.13 

'^•11^-0.12 


U.OU-o 03 

71±1 

>1000 

Eridanus 2 

56.0878 

-43.5332 

31.5 

22.9 

380 

-6.6±0.1 

l-isio.lo 

1 15Q+0.14 
i.oo_o 09 

169l^®8° 

Qg 

81±6 

1113 

Horologium 1 

43.8820 

-54.1188 

28.4 

19.5 

79 

-3.4±0.1 


1 “tl +019 
1-'^1-0.14 

3011.3 

< 0.28 


0.35 

Pictoris 1 

70.9475 

-50.2830 

17.3 

20.3 

114 

-3.1±0.3 

0.70l°;?^ 

0.881“:?^ 

2911? 

0 47+0-12 
-0.29 

78±23 

1.41 

Phoenix 2 

354.9975 

-54.4060 

13.9 

19.6 

83 

-2.8±0.2 

0.821°;?^ 

1091“:?® 

2611.? 

0 47+0-08 
' -0.29 

164T54 

1.81 

Indus 1 

317.2044 

-51.1656 

13.7 

20.0 

100 

-3.5±0.2 

0.84l°;5^ 

l.ZO_o 27 

qy+13.1 
-8.0 

< 0.38 


0.46 

Grus 

344.1765 

-50.1633 

10.1 

20.4 

120 

-3.4±0.3 

-1 qq + 0.74 
1.00_Q 26 

^^,^+0.85 
-L- < < -0.39 

^q + 29 . 8 
^^-13.6 

0 4]^+0-20 
'^•^-'--0.28 

4±60 

1.01 

Eridanus 3 

35.6897 

-52.2837 

10.1 

19.7 

87 

-2.0±0.3 

0-38lg:“ 

0 54+0-50 

14+12-5 
—2.6 

n 97+0.30 

< _o.l6 

83±36 

0.81 

Tucana 2 

342.9796 

-58.5689 

8.3 

18.8 

57 

-3.8±0.1 


9.83tj®® 

-L00_18.5 

0 39“*"°'^° 

u.aU-0 20 

107±18 

2.14 


Note. — Morphological parameters of the satellites are maximum a posteriori estimates with 68% (Icr) credible intervals, or limits 
Significance of detection in a (t-value) 

^ The uncertainty in distance modulus is estimated to be 0.1 —0.2 mag 
Half-light radius takes into account the ellipticity of the object via \/l — e multiplier 
The error on distance was not propagated to the physical size 
® Bayes factor for the elliptical vs circular model 

^ As this object is located very close to the CCD chip gap, its morphological properties should be treated with caution 


PSF magnitudes of star-like objects are given by 
the MAG_PSF output of SExtractor. As an indica¬ 
tor of star-galaxy separation we use the SPREAD_M0DEL 
parameter provided by SExtractor. This is a met¬ 
ric similar to psfmag-modelmag used by SDSS (see 
Fig. E}. A sensible selection thresho ld for bright stars 
would be [SPREAD MDDEL I < 0.003 (jPesai et al.l 120121: 
lAnnunziatella et al.ll2013l ). however for faint magnitudes 
this cut causes significant incompleteness in stars. There¬ 
fore, instead we choose to require 0 

|SPREAD_M0DEL| < 0.003 + SPREADERR.MODEL (1) 

This particular cut ensures that the stellar complete¬ 
ness remains reasonably high at faint magnitudes, while 
the contamination is kept low at the same time. The 
behaviour of 0.003 + SPREADERR_M0DEL as a function 
of magnitude shown in Eigure [T] explains why a fixed 
SPREADERR_M0DEL threshold is suboptimal. To assess the 
levels of completeness and contamination induced by our 
stellar selection, we use a portion of the DES-covered 
area of sky overlapp ing with the CEHTLS Wide survey 
(|Hudelot et al.ll2012l ). This is a dataset of comparable 
depth, for which morphological object classifications are 
provided. Eigure [2] gives the resulting performance of the 
stellar selection procedure in which Equation [T]is applied 
to both g and r-band catalogues. In particular, the Eig¬ 
ure gives completeness (black solid histogram) calculated 
as the fraction of objects classified as stars by CEHTLS 
(their CLASS_STAR>0.5) which are also classified as stars 
by our cuts applied to the DES data. Similarly, con¬ 
tamination can be gleaned from the fraction of objects 
classified as galaxies by the CEHTLS but as stars by our 
DES cuts (red dashed line). It is reassuring to observe 
low levels of contamination all the way to the very mag¬ 
nitude limit of the DES survey. At the same time, com¬ 
pleteness is high across a wide range of magnitudes and 
only drops to ~ 60% for objects fainter than r ^ 22. It 
is also worth noting that the star-galaxy separation cri¬ 
teria employed in this work may not be ideally suited for 
other studies, as they may have different requirements in 

^ http://I.usa.gov/lzHCdrq 


terms of the balance between the completeness and the 
contamination. 

In the stellar catalogues built using the procedure 
described above, the magnitudes are equivalent to the 
SDSS gri. Consequently, the extinction coefficients 
used are those suitable for the SDSS photometric sys- 
tem, while the dust reddening maps employed are from 
iSchlegel et al.l (|1998f ). Note that the depth of the result¬ 
ing catalogues varies somewhat across the DES footprint, 
but could be approximately estimated from the source 
number counts in g^ r, i filters. These number counts 
peak at magnitudes 23.7, 23.6, 22.9 in r, i correspond¬ 
ingly, indicating that the catalogues start to be signifi¬ 
cantly affected by incompleteness at somewhat brigher 
magnitudes g^23.5, r^23.4, i^22.7. 

To illustrate the quality of the resulting catalogue, Eig¬ 
ure [3] displays the density of the Main Sequence Turn- 
Off (MSTO; 0.2<g—r<0.6) stars on the sky. The den¬ 
sity of stars with 19<r<21 (corresponding to distances of 
10—25kpc) is shown in the green channel, more distant 
stars with 21<r<22.75 (corresponding to distances of 
25—56kpc) are used for the red channel, and the nearby 
stars with 17<r<19 (distances of 4—lOkpc) in the blue 
channel. This map is an analog of the ”Eield of Streams” 
picture bv iBeloknrov et ^ (|2QQ6l ). The density distribu¬ 
tion is very uniform thus confirming the high precision 
and the stability of the photometry as well as the robust¬ 
ness of the star-galaxy separation across the survey area. 
The map also reveals some of the most obvious overden¬ 
sities discovered in this work, at least two of which are 
visible as bright pixels in the Eigure. 

3. SEARCH FOR STELLAR OVER-DENSITIES 

To uncover the locations of possible satellites lurking 
in the DES data, we follow the approac h described in 
iKoposov et al.l (|2QQ8l ): IWalsh et aTT(j2QQ9[ ). In short, the 
satellite detection relies on applying a matched filter to 
the on-sky distributions of stars selected to correspond 
to a single stellar population at a chosen distance. The 
matched filter is simply a difference of 2D Gaussians, 
the broader one estimating the local background density, 
while the narrow one yielding the amplitude of the den- 
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Figure 4. Reticulum 2 satellite. Top row, panel: Density of all stars centered on the satellite’s location, dark shades of grey correspond 
to density enhancement. Red circles show the size of the apertures used to create CMDs given in the 1st and the 3rd bottom panels. 
Note the prominent over-density in the centre. Top row, 2'^^ panel: Density of stars selected using the isochrone mask shown in the 1st 
bottom panel. The background/foreground density variations have now subsided further as compared to the left-most panel. Top row, 
3'^^ panel: Density of galaxies in the vicinity of the satellite. No obvious over-density coinci dent with the object is visible. Top row, 4^^ 
panel: Distribution of Galactic dust reddening E(B — V) according to the measurement of ISchlegel et al.1 (|1998lb The extinction values 
corresponding to black and white colors in the panel are given in the panel title. Bottom row, panel: Colour-Magnitude distribution 
of stars inside the smaller of the two apertures shown in the 1st Top panel, centered on the object. Note the unmistakable overdensity 
corresponding to the Main Sequence and Main Sequence Turn-Off stars. The Red Giant Branch stars are also visible, albeit more scarce. 
There appears to be at least one possible Blue Horizontal Branch star member, as indicated by the good match between with isochrone 
(red dashed line) and the M92 BHB ridge-line (blue solid line). Bottom row, 2^^ panel: Comparison CMD built from random nearby 
stars outside the larger aperture, but covering the same on-sky area as in leftmost bottom panel. The comparison shows various stellar 
pile-ups, most prominently the disk dwarfs and objects near the detection limits, but features similar to that of the satellite’s populations 
are not observed. Bottom row, 3^^ panel: Satellite’s CMD for stars within the larger of the two apertures shown in the 1st top panel. At 
least 3 likely BHB members are visible. Bottom row, 4^^ panel: Hess difference for stars within and outside the aperture of r = 18^ Note 
the satellite’s MS, MSTO and RGB populations standing out clearly. Bottom row, 5^^ panel: Density of nearby background stars in the 
colour-magnitude space (Hess diagram). 


sity peak at the location of the satellite. 

We start by taking a catalogue of sources classified as 
stars. A sub-set of these is then carved out with either a 
set of colour-magnitude cuts or with an isochrone mask 
offset to a trial distance modulus. Then a 2D on-sky 
density map of the selected stars is constructed, keeping 
the spatial pixel sufficiently small, e.g. V on a side. At 
the next step, the density map is convolved with a set 
of matched filters (described above) with different inner 
and outer kernels. Finally, these convolved maps are 
converted into Gaussian significance maps and the most 
significant over-densities are extracted. 

In the analysis presented here, we have used the mask 
based on the PARSEC isochrone (jBressan et al.l l2Q12f ) 
with an age of 12 Gyr and metallicity of [Fe/H] = —2 0 , 
spread out by the characteristic DES photometric error 
as a function of colour and magnitude. The trial distance 
moduli explored in the search range from m—M=15 to 

® http://stev.oapd.inaf.it/cmd 


m—M=24. The inner kernel size is allowed to vary from 
V to 10'. After extracting the most significant over¬ 
densities all detections are cross-matched with the list 
of pos itions of known dwarf ga laxies from iMcConnachi^ 
(j2Q12[ ). nearby LEDA galaxies (jPatnrel et al.ll2QQl ) with 
radial velocities l ess than 2000 km/s, globular clusters in 
the catalogue of iHarrisI (jl996f ) as well as the globular 
clusters listed in Simbad. 

The list of over-densities produced is then cleaned for 
duplicates, e.g. over-densities detected with more than 
one distance modulus and/or kernel size. We have also 
removed detections near the very edge of the survey’s 
footprint. The remaining candidates are ranked by their 
significance and eye-balled. Among the 9 objects pre¬ 
sented in this paper, 5 are at the very top of the ranked 
list with significances well above lOcr, while a further 4 
are somewhat lower in the list (but with significances 
of order of lOcr). We firmly believe that all objects pre¬ 
sented below are genuine new satellites of the Milky Way. 
Apart from the over-density significance, our inference 
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Figure 5. Eridanus 2. See Figured for detailed description of the panels. This is the most distant of the detected satellite. Note the 
faint over-density of objects classified as galaxies at the location of the satellite. These are mis-classified blended stars unresolved by 
DES+SExtractor. The CMD shows prominent RGB and BHB with a handful of possible young Blue Loop stars. 


is based on a combination of additional factors, includ¬ 
ing their morphological and stellar population properties. 
We describe this overwhelming observational evidence in 
detail in the next Section. The list of objects and their 
properties is given in Table [T] and discussed in detail in 
the next Section. 

4. PROPERTIES OF THE DETECTED OBJECTS 

While our search for stellar over-densities returned a 
large number of promising candidates, in this Paper we 
concentrate on a particular gold-plated sample of objects 
whose nature can be unambiguously established with the 
DES data alone. Let us briefly review the sanity checks 
carried out to corroborate the classification. First, each 
of the nine objects possesses a significant stellar over¬ 
density, in fact, the significance is in excess of lOcr in 
eight out of nine cases. Second, no obvious background 
galaxy over density is visible, thus reducing dramatically 
the possibility that the objects are spurious density peaks 
caused by mis-classified faint galaxies. Third, the dust 
maps in the vicinity of each of the satellite do not show 
any strong features that can be linked to the stellar over¬ 
density. Most importantly, distributions of stars in the 
colour-magnitude plane reveal sequences corresponding 
to coeval populations at the same distance. In all cases, 
not just one particular stellar population is visible but 
several, for example, MSTO+RGB+BHB. 

Figures |3][T2] give the visual summary of the properties 
of each of the detected satellites and serve to substanti¬ 
ate our classification. Each of these diagnostic plots com¬ 
prises of 9 panels: 4 in the top row are concerned with 
the on-sky distribution of the object’s stars, while the 5 
in the bottom row give details of its stellar populations. 


In particular, the density of all objects classified as stars 
is shown in the top left panel. The second panel in the 
row gives the density of all stars selected with the best- 
fit isochrone mask (shown in the left bottom panel). For 
comparison, the density of objects classified as galaxies is 
presented in the third panel, an d the galactic extinction 
map from iSchlegel et al.l (|1998[ ) in the fourth. 

The first panel in the bottom row displays the colour- 
magnitude distribution of stars selected to lie in the small 
aperture (whose size is indicated by the small red circle 
in the bottom right corner of the top left panel) cen¬ 
tered on the object. Red solid lines are the bound¬ 
aries of the isochrone mask used to select the proba¬ 
ble satellite member stars, while the isochrone itself is 
shown as a red dashed line. The horizontal branch (HB) 
in the model isochrone does not extend to sufficiently 
blue g — r colours, therefore we complement the theoret¬ 
ical HB with a fiducial HB (blue solid line) built using 
data for the metal-p oor old globular cluster M92 from 
[Bernard et al.l (|2Q14D . The object’s CMD can be con¬ 
trasted to the foreground/background CMD shown in the 
second bottom panel. This is built with random nearby 
stars located within large annulus (typical inner radius 
0.2° and outer radius of 1°) and sampling the same num¬ 
ber of stars as that in the left panel. The aperture chosen 
to create the CMD in the left panel is sufficiently small 
in order to minimize the contamination by non-member 
stars. This, however, implies that some of the rarer stel¬ 
lar populations (such as BHBs for example) are not fully 
sampled. A more complete view (albeit with a higher 
contamination) of the satellite’s stellar populations can 
be seen in the third bottom panel for which a larger aper- 
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ture is used (large red circle in the top left panel). In 
most cases, this particular view of the CMD strength¬ 
ens the detection of the BHB population. The last two 
panels in the bottom row are Hess diagrams, i.e. density 
distributions in the colour-magnitude space. In particu¬ 
lar, the fourth panel is the Hess difference between the 
stars inside a small aperture centered on the satellite (the 
aperture size is shown in the panel’s title) and a large 
area well outside the satellite’s extent. Finally, for com¬ 
parison, the Hess diagram of the background/foreground 
population is shown the fifth panel of the bottom row. 

All objects we report in Tabled] have passed the sanity 
checks described above with flying colours. Moreover, 
some of the new satellites are actually visible in the DES 
image cutouts as described in the sub-section 14.21 

The naming convention we have adopted uses the con¬ 
stellation of the object together with the corresponding 
arabic numeral. This differs from the customary bimodal 
naming convention for globular clusters and dwarf galax¬ 
ies, that leads to the name dependance on the knowledge 
of whether the object is a cluster or a galaxy. 

4.1. Satellite modelling 

Having identified the objects of interest, we proceed 
to measuring their distances, luminosities and structural 
parameters. First, the distances are determined using 
the available BHB stars if present. We have used both 
the theoretical HB isochrone corresponding to the old 
metal-poor population mask used for the selection of the 
likely satellite memb ers and the M92 HB fiducial from 
[Bernard et al.l (j2Q14f ) converted into the SPSS photomet - 
ric system using the equations from iTonrv et al.l (|2Q12f ). 
The distance moduli implied by these HBs are consis¬ 
tent within 0.1 mag. To verify the distance accuracy we 
also fitted the Hess diagra m of each object with a stellar 
population model (see e.g. iKoposov et al.l[2QT(i for more 
details). For all objects (excluding Tucana 2) this gave 
us similar distance modulus differing by less than 0.1 - 
0.2 mag. For Tucana 2 we use the distance implied by 
the Hess diagram fit. 

To measure the morphological properties of all the de¬ 
tected objects, we have modelled the on-s ky 2D distri¬ 
butio n of the likely member stars (see e.g. [Martin et al.l 
[2008L for similar approach). For each object, only stars 
located within the CMD mask are used for the mor¬ 
phological analysis. The model used to describe the 
density distribution of stars is a mixture of the back¬ 
ground/foreground stellar density (assumed to be con¬ 
stant across a few degrees around the object) and a 
rotated elongated exponential model for the satellite’s 
stars. Therefore, the probability of observing a star at a 
particular position on the sky with coordinates a and 5 
is: 


P(a, ^1^) = (!-/) P6^(a, d\0) + fPobj{a, d\0) 


Pbg{a,5\9) 


1 

D cos^c 


Pohj{a,5\0) 


27r(l 


1 _/ Dell \ 

e) COsJc<^^aj|- \ (^maj J 


and 

Dell 


cos(PA) sin(PA) \ 
(1-e) ^ (1-e) ^ 

— sin(PA) cos(PA )J 



- ac) cosSc\ 

s-s, ) 


Here Q is the area of the analysed part of the data and 
^ is a shorthand notation for all model parameters, of 
which there are 6 in total. These are: / - the fraction of 
objects belonging to the object, Omaj ~ the major axis 
(the exponential scale length), e - ellipticity of the object, 
PA positional angle of the major axis, and finally, ac Sc 
- the position of the center of the object. 

The above model is fit to the data assuming uninfor¬ 
mative priors on ac^Sc, hat priors on e and PA, (i.e. e ^ 
Uniform(0, 0.8), PA ^ Uniform(0, 180)) and the Jeffreys 
priors on / and Omaj^ namely / ^ Beta(0.5, 0.5) and 
P(amaj) oc ^ (see e.g. iMackavl 120021: [Gelman et al.l 

120141 ). The resulting posterior is sampled using the Affine 
invariant Ensemble sampler (jGoodman fc Wear3 [20101: 
[Foreman- Mackey et al.|[2M^ . The values of parameters 
and the associated uncertainties are given in Tabled] To 
check whether the ellipticity of an object is significant, 
the B ayes Factor is calculated using the Savage-D ickey 
ratio (jVerdinelli fc Wassermannl [19951 : iTrotta] 120071 ). As 
the Tabled] reports (final column), only Reticulum 2 and 
Eridanus 2 show strong, statistically significant elonga¬ 
tion. Other objects are consistent with being circular. 

Having modelled the morphology of each satellite, the 
posterior distributions of their structural paramereters 
are used to estimate the total number of stars belonging 
to the satellite. The satellite luminosities are calculated 
only with stars with r<23 falling within the isochrone 
mask (we have used a brighter magnitude limit to avoid 
incompleteness due to the star/galaxy misclassification). 
Then, the number of stars inside the isochrone mask is 
converted into the stellar mass and the corresponding to¬ 
tal luminosity assuming a n old metal-poor isochrone and 
the KTG/Chabrier IMF (jKroupa et al.l [19^ iChabriei] 

inni). 

Figure dS] is a gallery of ID azimuthally averaged stellar 
density profiles for each of the satellites. Red lines are the 
best fit exponential models. Before we discuss the size- 
luminosity distribution of our sample in sub-section 14.3[ 
let us have a closer look at the individual properties of 
each of the satellites. 


4.2. Notes on individual ohjeets 

• Reticulum 2 (see Figure [4|) - With the significance 
in excess of 40, this is the most obvious previously 
unknown object in the DES field of view. The 
CMD shows clear MS and MSTO together with be¬ 
ginnings of a RGB. Reticulum 2 is clearly very elon¬ 
gated: the measured axis ratio is around 0.6. It is 
one of the two objects in the sample presented here, 
for which the flattening of the stellar density is sta¬ 
tistically significant. Its luminosity is My = — 2.7 
and the half-light radius of 30 pc, which brings 
it just outside the cloud of classical globular clus¬ 
ters into either extented globular clusters or dwarf 
galaxies. It is impossible to classify this satellite 
securely with the data in hand, therefore we shall 
designate this and similar objects as ’’ultra-faint 
satellites”. Clearly, kinematic and chemical analy¬ 
sis is required before the nature of this object can 
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Figure 6. Horologium 1 satellite. See Figure [4] for detailed description of the panels. The CMD reveals an obvious MSTO and RGB. 
There are also several likely BHB members. 
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Figure 7. Pictoris 1 satellite. See Figure [4] for detailed description of the panels. Pictoris 1 resembles Horologium 1 (see Figure [6] albeit 
located at a slightly larger distance. Therefore the MSTO is barely detected as it kicks in right around the DES limiting magnitude. 
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Figure 8. Phoenix 2 satellite. See Figure [4] for detailed description of the panels. The satellite looks similar to Horologium 1 and Pictoris 
1, but with a noticeable Blue Straggler sequence. 
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Figure 9. Indus 1 satellite. See Figure [4] for detailed description of the panels. Located at a similar distance to Horologium 1, Pictoris 1 
and Phoenix 2, Indus 1 is more co mpa ct and boasts only one possible BHB star. Curiously the satellite is visible in the DES image-cutout 
presented in the companion Figure fl^ 
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be established. It is reasonable to associate the low 
axis ratio measured for this object with the effects 
of the Galactic tides. However, at the moment, we 
do not have any hrm evidence of extra-tidal ma¬ 
terial outside the satellite. The distance estimate 
for Reticulum 2 is 30 kpc as judged by a couple of 
BHB stars as well as the prominent MS and RGB. 
Note that eight other satellites are all substantially 
further away. 

• Eridanus 2 (see Figure [5j). This is the second most 
signihcant and the most distant object in our list. 
With its location on the periphery of the Galaxy 
at ^ 380 kpc, luminosity of My = — 6.6 and the 
size of ^ 170 pc, this appears to be a twin of the 
Leo T dwarf galaxy (jirwin et al.ll2007l ). Similarly 
to Leo T, Eridanus 2 is clearly visible in the DES 
color images (Figure [15]) as a group of faint blue 
stars embedded in a blue low surface-brightness 
cloud. In contrast to Leo T though, Eridanus 2 
shows signihcant elongation with a measured ellip- 
ticity of ^ 0.4. Its color-magnitude diagram boasts 
a strong RGB and a very prominent red clump. 
Additionally, there appears to be a group of bright 
blue stars which could possibly be interpreted as 
young blue-loop stars (see Fignre [T6l) . similar to 
those found in Leo T (jWeisz et al.ll201^ . The HB 
stars of Eridanus 2 are just at the magnitude limit 
achievable with the single-frame DES photometry. 
The possible presence of young stars and large dis¬ 
tance to the Eridanus 2 galaxy suggest that the 
satellite is a viable candidate for follow-up HI ob¬ 
servations. The color images of the galaxy (see 
both panels of Figure [15] also show an curious fuzzy 
object which can be interpreted as a very faint glob¬ 
ular cluster. The object is almost at the center of 
the galaxy (a, (5) = (56.0930,—43.5327) and is very 
diffuse with a size consistent with a few parsecs as¬ 
suming it is at distance of Eridanus 20. If any of 
the faint diffuse objects superposed with Eridanus 2 
are indeed globular clusters, this would make this 
satellite the faintest galaxy to possess stellar sub¬ 
systems! 

• Horologium 1 (see Figure [6]). This is the third most 
obvious object within the DES held-of-view with a 
signihcance of ^ 30. Horologium 1 is one of 7 satel¬ 
lites all lying at approximately the same distance 
of ^100 kpc. Both its Hess diagram and the GMD 
clearly show a prominent red giant branch as well 
as the group of the BHB stars. With a luminosity 
of My = —3.5 and a size of 30 pc it sits within the 
uncertain territory where extended star cluster and 
faint dwarf galaxies start to overlap. We class it as 
an ’’ultra-faint satellite”. 

• Pictoris 1 (see Figure [3. Another high signihcance 
ultra-faint satellite located at ^ 115 kpc and pro¬ 
jected near the LMC on the sky. The GMD shows 
a clear RGB and a few stars on the HB. Its size is 
^ 30 pc and its luminosity is My = —3.2. 

^ There is another compact but high surface brightness round 
object looking consistent with a globular cluster (a, (5)= (56.0475,- 
53.5213) but it is probably a background galaxy as even its out¬ 
skirts do not resolve into stars 


• Phoenix 2 (see Figure IH|. An ultra-faint satellite 
at ^ 80 kpc with the half-light radius of 26 pc and 
luminosity of My = —2.8. Apart from a clear 
MSTO and RGB, one or two BHB stars, there is 
potentially a large number of BS stars extending 
to g—r<0, r^22. 

• Indus 1 (see Figure [9]). Similar to the bulk of our 
sample, this ultra-faint satellite is at 100 kpc with 
a half-light radius of 36 pc and the luminosity of 
My = —3.5. Interestingly, Indus 1 is conspicuous 
in the DES false-color image (see Figure [T4|) . 

At the time of the submission of this manuscript, 
the discovery of t he same satellite was announced in 
iKim et ^ (j2015[ ) where the object’s size and lumi¬ 
nosity are quoted to be signihcantly lower than the 
values we obtain. Usi ng imaging data m uch deeper 
than available to us, iKim et akl (j2015[ ) report the 
detection of a pronounced mass segregation, which 
(if true) would argue strongly in favour of the glob¬ 
ular cluster nature of the object. Moreover, the dis¬ 
agreement in the inferred total luminosity can be 
explained by the difference in depth between the 
two datasets in the presence of a depletion at the 
low-mass end of the mass function. With regards to 
the extent of the object, we note that it is possible 
that our size determination has been biased high by 
blending in the central parts of the object . How¬ 
ever, the small size reported by IKim et aTI (j2Q15[ ) 
could, at least in part, be the consequence of the 
choice not to ht the stellar density at lager radius, 
as it is classihed as “extra-tidal” by the authors. 

• Grus 1 (see Figure [TOj). The satellite is at 120 kpc, 
has a luminosity of —3.4 and a size of 60 pc. Given 
its luminosity and size, the satellite should per¬ 
haps be classihed as a dwarf galaxy. However, care 
should be taken when interpreting the measured 
structural parameters, as the object sits right at 
the CCD chip gap. 

• Eridanus 3 (see Figure El). The satellite is at 
90 kpc with the smallest half-light radius in the 
sample, of only 13 pc and the luminosity of My = 
—2.0. Eridanus 3 is also easily noticeable in the 
DES image (right panel of Figure [T4|) . 

• Tucana 2 (see Figure [T^ . This satellite can per¬ 
haps be classihed as an ultra-faint dwarf galaxy 
with less ambiguity than many others in the sam¬ 
ple, due to its more substantial luminosity of My = 
—4.4 and the large size of 160 pc. The extent of the 
object is immediately obvious from the map of se¬ 
lected stars (second panel in the top row). The 
GMD demonstrates a strong MSTO, a sub-giant 
branch. A few blue objects with colors consis¬ 
tent with the BHB stars are probably foreground 
contaminants as their magnitude contradicts the 
distance modulus implied by the subgiant branch. 
The estimated distance is ^ 60 kpc 

4.3. Size-luminosity relation 

Figure [17] shows the locations of the discovered satel¬ 
lites on the plane of half-light radius (1.67a^ajVl — e) 
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Figure 10. Grus 1. See Figure [4]for detailed description of the panels. Grus 1 is situated at a distance similar to that of Horologium 1, 
Pictoris 1, Phoenix 2 and Indus 1. The central stellar overdensity is not as prominent as in the previously discussed objects. This is partly 
explained by the fact that Grus 1 falls right into the GGD chip gap. Nonetheless, the GMD shows off such familiar features as an RGB 
and a well-populated HB. 


and intrinsic luminosity {My). Together with the ob¬ 
jects in our sample (red filled circles), several other stellar 
systems are displayed. These include globulars clusters 
(GCs, black dots), MW dwarf galaxies (black unfilled cir¬ 
cles), M31 dwarf galaxies (black empty triangles), as well 
as extended objects wi th half-light radius le ss than 100 pc 
from the catalogue of iBrodie et al.l (|2011[ ). Let us start 
by pointing out that all groups of objects in this compi¬ 
lation are affected by selection effects. Most importantly, 
the detectablility of all objects is dependent on their sur¬ 
face brightness. For objects detectable as resolved stellar 
systems in surveys similar to the SDSS, the parameter 
space b elow the surface-brig htness limit of /i ~ 31 mag 
deg“^ (jKoposov et al.l I^OOSl ) is unexplored (see dashed 
lines of constant surface-brightness in the Figure). Ad¬ 
ditionally, there are constraints set by blending for star 
clusters, which deteriorate with increasing distance. 

Notwithstanding various selection effects, star clusters 
and dwarf galaxies seem to occupy distinct regions of the 
size-luminosity plane. Brighter than My = — 6 , star clus¬ 
ters seem to rarely exceed half-light radius Vh ^ 60 pc, 
while no dwarf galaxy is smaller than 100 pc. In other 
words, there seems to be a size gap between star clusters 
and dwarf galaxies. At fainter magnitudes, the picture 
is a lot more muddled. Not only are there objects firmly 
classified as dwarf galaxies with size under 100 pc (for ex¬ 
ample Pisces II and Canes Venatici II), but also several 
ultra-faint satellites, namely Segue 1, Segue 2 and Bootes 
2 have sizes and luminosities comparable to faint ex¬ 
tended globular clusters such as Palomar 4 and Palomar 
14. What sets Segue 1,2 and Bootes 2 apart from sim¬ 
ilar looking star clusters are their chemical abundances. 


All t hree objects appear to have large metallicity spreads 
(e.g. iFrebel et al.l |2Q14 iKirbv et al.l 120131: iKoch fc Richl 
l2014f ). previously unnoticed in star clusters. 

Guided by the approximate shapes of the distribu¬ 
tion of star clusters and dwarf galaxies in size-luminosity 
plane, we can conclude the following. Using the mea¬ 
sured morphological parameters only, out of the 9 satel¬ 
lites, only three, namely Eridanus 2, Tucana 2 and Grus 
1 can be classified as dwarf galaxies unambiguously, 
without having to resort to follow-up spectroscopy. We 
classify the remaining six, i.e. Eridanus 3, Phoenix 2, 
Horologium 1, Pictoris 1, Reticulum 2 and Indus 1 sim¬ 
ply as ultra-faint satellites with the hope that their true 
nature is revealed through detailed abundance studies in 
the near future. 

It is interesting that in the size luminosity plane the 
ambigiously classified objects in our sample seem to join 
the sequence of objects observed in other galaxies la- 
bell ed as faint-fnzzies ( EE) or Extended Clusters (EG) 
(see Larsen et al.]l2QQll: IBrodie et al.ll2Qlll: lEorbes et al.l 
|2Q13[ ). These objects have sizes of tens of parsecs and 
luminosities of My ^ — 5. While our objects are signifi¬ 
cantly fainter, one should take into account that most 
HST surveys of other galaxies are incomplete for ob¬ 
jects with My < — 5. The possible connection between 
EE/EC objects and the ultra-faint satellites is particu¬ 
larly intriguing given the likely association between some 
of our objects and the Magellanic Clouds (see Section [5]). 

Strikingly, while most of the new satellites announced 
here are rather faint, none of them lies particularly close 
to (with the exception of Tucana 2), or indeed below, 
the nominal surface brightness limit of /i ^ 31. This 
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Figure 11. Eridanus 3. See Figure[4]for detailed description of the panels. A very compact satellite that falls in the distance range shared 
by Horologium 1, Pictoris 1, Phoenix 2, Indus 1 and Grus 1. The satellite’s CMD sequence is rather sparse. There are hints of a turn-off, 
BS stars and an RGB. Despite its less convincing GMD, Eridanus 3 is a genuine satellite as can be seen from the DES image cut-out shown 
in Figure [T4| 


means that these satellites are not very different from 
those discovered by Sloan (only somewhat more distant 
at fixed luminosity), and, therefore, significant fraction, 
but not all of them (as the ob j ects lie very close to th e lu- 
minosity vs distanc e limits of iKoposov et alJ (|2QQ8[) and 
I Walsh et al.l (|2QQ9[ )) would have very likely been identi¬ 
fied by the SDSS and the VST ATLAS had they fallen 
within the footprints of these surveys. The relatively 
high surface brightness values of the DES satellites must 
imply that we have not tapped into a large supply of 
nearly invisible dwarfs predicted by some of the semi- 
analytic models. Instead, the discovery of a large num¬ 
ber of satellites from a small area as compared to the 
previously observed sky suggests a peculiar spatial dis¬ 
tribution of these objects. This is the focus of the next 
Section. 

5. SPATIAL DISTRIBUTION 

Let us start by inspecting the overall distance distri¬ 
bution of the DES satellites. Eigure [18] gives the ob¬ 
ject’s luminosity as a function of Galacto-centric radius 
RgC‘ It is immediately obvious that i) our satellites are 
some of the faintest of the distant MW satellites known, 
and ii) a large faction of them cluster around the dis¬ 
tance Rgc ^ 100 kpc. In other words, the satellites as a 
group lie not far from the LMC {Rgc ^ 50 kpc) and just 
behind the orbit of the Small Magellanic Cloud (SMC, 
Rgc ^ 70kpc). 

A more direct test of a possible association with the 
Magellanic Clouds is shown in Eigure [191 Here we plot 
the distance to the LMC against the distance to the 
Galactic centre for the DES satellites as well the known 


MW dw yf galaxies (iMcConnachid 120121) and globular 
clusters (|Erancis fc Anderson! 20141) . In this space, the 
objects that lie above the one-to-one diagonal line of the 
plot are closer to the LMC than to the Galactic centre. 
Note that all our objects lie above this demarcation line. 
While surprising at first, this is in fact expected due to 
the particular sky coverage of the DES survey (see Eig¬ 
ure However, even if the bunching around the Mag¬ 
ellanic Clouds is set by the footprint, the total number 
of the satellites uncovered in such small area is alarm¬ 
ing, especially when compared to the Southern portions 
of the SDSS, or VST ATLAS. 

To test whether there is a statistically significant over¬ 
density of satellites in the vicinity of the LMC we carry 
out the following calculation. The 3D position vector 
of each MW satellite, including the newly discovered 
ones, but excluding the LMC itself, is rotated randomly 
50000 times with respect to the Galactic center (GC) 
while keeping fixed its Galacto-centric radius. After each 
rotation, the distance to the LMC Rlmc is computed 
and stored. The blue/white background density in Eig¬ 
ure [19] shows the average expected density of objects in 
the (Rlmc, Rgc) space from the random 3D position 
reshuffles described above. Under the assumption of an 
isotropic distribution of satellites in the Galaxy, our pro¬ 
cedure would yield a reasonable estimate of the expected 
number of satellites in each portion of the (Rlmc, Rgc) 
plane. Darker shades of blue show areas where more 
objects are expected, while pale blue and white pixels 
correspond to less populated corners of the parameter 
space. 
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Figure 12. Tucana2. See Figure [4]for detailed description of the panels. Tucana 2 is the closest of the group of 7 satellites at around 60 
kpc. It is considerably fluffier in appearance which explains its lower signiflcance. The strongest feature in the CMD is the MSTO. 


According to the simulation described above, the entire 
group of the satellites announced in this paper (excluding 
the distant Eridanus 2) is located in the under-populated 
area of the {Rlmc^ Rgc) space, as indicated by the 
pale blue background density. This can be interpreted 
as a detection of an over-density of satellites around 
the LMC. To quantify the significance of this detec¬ 
tion we perform the following frequentist hypothesis test. 
We carve the LMC-dominated area in the (Rlmc^Rgc) 
space: Rlmc < Rgc and Rlmc < 100 kpc. Then the to¬ 
tal number of satellites in the LMC-dominated area (ex¬ 
cluding LMC itself) is computed. The actual observed 
number {N = 14) is compared to the number of objects 
in the marked region as given by the random reshuffles of 
the 3D position vectors. The distribution of the possible 
number of satellites in the marked area is shown in the 
inset of the Figure, and the observed number is given by 
the red vertical line. The distribution obtained assuming 
isotropy is narrow and peaks around = 10, indicating 
the presence of a satellite overdensity in the DES data. 
However, the formal significance level of this result is not 
very high at ^ 94/^ while the estimated excess is also 
quite modest, i.e. of order of 4 objects. 

Having established a possible connection between the 
Magellanic Clouds and the DES satellites, we study 
the distribution of the new objects in the plane of the 
LMC’s orbit. Figure [20] shows the MW satellites (GCs, 
dwarfs and DES satellites) projected onto the plane de¬ 
fined by the vector of the LMC’s velocity as given in 


® The boundary in the Rlmc vs Rqc space used for the test is 
chosen to be as simple as possible to avoid artiflcially inflating the 
p-value. 


iKallivavalil et ah! (j2013[ lFI. We also show the forward- 
and backward-integrated orbits of the LMC in two MW 
potentials with different DM halo concentrations (c = 10 
and c = 25). Coloured regions correspond to the frac¬ 
tion of the Galactic volume covered by the SDSS, VST 
ATLAS and DES surveys as projected onto the LMC’s 
orbital plane. The three slightly over-lapping surveys 
have together covered a large portion of the sky. How¬ 
ever, some lacunae still remain, most notably directly in 
front of the LMC-SMC pair as judged by the vector of the 
LMC motion (black arrow). As expected, the DES satel¬ 
lites bunch tightly outside the SMC’s orbit. However, 
while the objects have similar in-plane coordinates, not 
all of our satellites actually lie close to the orbital plane 
of the LMC. The distance from the LMC’s orbital plane 
is shown in Figure |2T1 Only, Reticulum 2, Horologium 
1 and Eridanus 3 have small heights above the plane 
\zlmc\ < 10 kpc. The distribution of the other 6 objects 
do not show any strong alignment with the LMC’s orbit. 

The accretion of the Magellanic system has left be¬ 
hind a trail of neutral hydrogen. Figure [22] compares 
the positions of the DES satellites to the distribution of 
HI in the Magellanic Stream as traced by iPutman et ah! 
(j2QQ3[) . Curiously, the DES satellites seem to avoid the 
high column density portions of the stream. The spa¬ 
tially distinct behavior of the gas and the stars is simi¬ 
lar to that seen around M31 (jLewis et al.ll2Ql^ . There 

® The transformation from the Galactocentric coordi¬ 
nates (X,Y,Z) (X is positive towards the Galac tic antic en- 
ter) to Xlmc^YlmCi Zlmc shown in Figures [20] and [21] 
is deflned by Xlmc = 0.10490391 X - 0.99448236 T + OZ; 
Wmc = 0.14676983 X + 0.015482154 Y + 0.98904950 Z; Zlmc = 
-0.98359229 X - 0.10375516 Y + 0.14758415 Z. 
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Figure 13. Azimuthally averaged ID density profiles of the detected satellites (black curves) with the associated uncertainties. Red curves 
represent the best-fit exponential density profiles. Note that the actual fits has been done using unbinned 2D data. 



Figure 14. DES false-colour composite images of the two compact satellites, Indus 1 and Eridanus 3. Each image is 4.3^x4.3h The R,G,B 
channels correspond to the z, r, g bands. 


is however possible overlap between Horologium 1, Eri¬ 
danus 3, Phoenix 2 and Tucana 2 with the high veloc¬ 
ity clouds associated with the Stream. Additionally, the 
particular projection used in the Figure emphasizes the 
fact that there appears to be two groups: one contain¬ 
ing Pictoris 1, Eridanus 2, Reticulum 2, Horologium 1 
and Eridanus 3, all lying on the EMC side of the Stream. 
Meanwhile, Phoenix 2, Tucana 2, Grus 1 and Indus 1 
can be found on the SMC side. Interestingly, the three 
objects that have the highest probability of being aligned 
with the orbital plane of the LMC (see Figure [21]), i.e. 


Reticulum 2, Horologium 1 and Eridanus 3 are the satel¬ 
lites closest to the Magellanic Stream. 

6. DISCUSSION AND CONCLUSIONS 

In this paper, we have presented the discovery of 9 new 
ultra-faint Milky Way satellites using the data from the 
Dark Energy Survey. Based on the morphological prop¬ 
erties, 3 of the satellites are definite dwarf galaxies, while 
the exact type of the other 6 is more ambiguous. The 
nearest object in the sample shows a significant elonga¬ 
tion indicative of the influence of the Galactic tides. The 
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Figure 15. DES false-colour composite image of the Eridanus 2 dwarf galaxy. Left: 13'xl3' field around the galaxy, which is visible 
as a clear concentration of blue stars as well as the low surface brightness emission due to unresolved stars. Right: The central regions 
(6.5'X6.5') of the Eridanus 2 dwarf galaxy. The very low surface brightness round spot is clearly visible near the center and is possibly a 
globular cluster associated with the galaxy. A brighter elongated blue object below is more likely to be a background source. 


CMD r<4.0’ 





— Age=12Gyr[Fe/H]=-2 

— Age=250Myr [Fe/H]=-0.75 



-0.5 0.0 0.5 1.0 1.5 2.0 
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Figure 16. Evidence of young stellar populations in the color- 
magnitude diagram of the Eridanus 2 dwarf galaxy. As well as the 
old metal-poor stars (12 Gyr, [Fe/H] = —2 isochrone, red curve), 
the dwarf galaxy seems to show a significant Red Clump and a 
handful of Blue Loop stars indicative of recent star-formation as 
revealed by the metal-rich 250 Myr isochrone (blue curve). 

most distant object discovered appears to be a dwarf 
galaxy located at the very edge of the Milky Way halo, 
around ^ 380 kpc. This remote dwarf shows signs of re¬ 
cent star formation and possibly even hosts a faint glob¬ 
ular cluster! The discovery of a large number of MW 
satellites in a small area around the LMC and the SMC 
is suggesting that at least a fraction of the new objects 
might have once been a part of the Magellanic group. 


The idea that some of the satellite galaxies of the Milky 
Way may be associated with the Magellanic Clouds is 
quite old. Let us re-visit the earlier work in the light of 
t he discoveries in thi s paper. 

iLvnden-Belll (|1976f) first suggested the idea of a Greater 
Magellanic Galaxy. His paper predates modern ideas 
on galaxy formation, and so his picture is rather differ¬ 
ent from today’s concensus. He envisaged dwarf galax¬ 
ies condensing from gas torn out of the Greater Magel¬ 
lanic Galaxy via a tidal encounter with the Milky Way 
Galaxy. He noted that Draco and Ursa Minor lie in 
the Magellanic Plane, almost opposite the Magellanic 
Glouds in the Galactocentric sky, and were candidates 
for debris from the break-up of the Greater Magellanic 
Gala xy. This speculatio n received some confirmation 
when ILvnden-Belll (|1982D showed that the elongation of 
Draco and Ursa Minor is along th e Magellanic Stream. 
iLvnden-Bell fc Lvnden-BdH (|1995f ) later suggested that 
Sculptor and Carina can also be added to the associa¬ 
tion of galaxies in the Magellanic Stream, though their 
membership is more tentative as their orientations are 
not along the Stream. Nonetheless, the evidence is at 
least suggestive for a linkage between Draco, Ursa Minor 
a nd the Magellanic Clouds . 

iD’Onghia fc Lak^ (|2QQ8f ) argued that the Magellanic 
Clouds were the largest members of a group of dwarf 
galaxies that fell into the Milky Way halo at a com¬ 
parati vely late time. This was inspired by iTnllv et ahl 
(|2QQ6f) observational discovery of a number of nearby as¬ 
sociations of dwarf galaxies just outside the Local Group 
at distances of ^ 5 Mpc. These groups have masses 
1 - 10 X 10^and typically contain ^ 5 dwarfs 
comparable to in mass to the larger dwarf spheroidals, 
as well as an unknown number of fainter satellites. The 
sizes of the dwarf associations are a few hundred kilopar- 
secs, and the velocity dispersions a few tens of kms“^. 
They are bound, but not in dynamical equilibrium as the 
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Figure 17. Luminosity as a function of half-light radius for star clusters and galaxies. M W globular c l usters (black dots, 
IFrancis &: AndersonI (12014111. extended obje cts with half-light radius less than 100 pc (grey dots , from IBrodie et al.l (|20111 11. MW dwarf 
satellites (black circles. [McConnachid (|201^ 11. M31 dwarf satellites (black triangles. IMcConnachid (|2012I1 1 and the DES satellites (red filled 
circles) are shown. Dashed grey boundaries are curves of constant surface brightness = 26 mag deg“^ and = 31 mag deg“^. Grey 
hatched region is the portion of the parameter space currently unaccessible to surveys like SDSS and DES. Note that the DES satellites do 
not lie too close to the low surface-brightness limit and, therefore, could have been detected by SDSS. Most of the DES objects (apart from 
Eri 2, Tuc 2 and Gru 1) lie in the region right between the extended globular clusters and faint dwarf galaxies. Thus, their classification 
remains ambiguous with the data in hand. 
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Figure 18. Luminosity as a function of Galactocentric distance 
for globular clusters (black dots), MW dwarf galaxies (black un¬ 
filled circles) and the DES satellites (red filled circled). Note a 
large group of the satellites congregating just outside the SMG 
orbit (vertical dashed line). 

crossi ng time is a sizeable fract ion 80%) of the Hubble 
time. iD’Onghia fc Lak^ (j2QQ8f 1 speculated that the LMC 
may have been the largest member of such a dwarf asso¬ 
ciation. When this Magellanic Group fell into the Local 
Group, it provided seven of the bright satellites of the 
Milky Way (LMG, SMG, Draco, Sex, Sgr, UMi, Leo II). 
They argued that subhalos accreted singly may remain 
dark, whilst subhalos accreted in groups may light up. 
This is because the low virial temperature in the parent 


halo allows gas to cool and be re-accreted onto subhalos, 
though this argument was not corroborated by hydro- 
dynamical simulations. Altho ugh it is unclear whether 
the seven dwarfs proposed by iD’Onghia fc Lak^ (j2QQ8l ) 
really all were part of a Magellanic Group, it is certainly 
true that gro up infall occnp frequently in numerica l sim¬ 
ulations (e.g. ILi fc Helmill200l IWetzel et ^ lM^ and 
S O the underl y ing m otivation is sound. 

iSales et al ] (l2fm]] studied the orbits of the analogues 
of the Magellanic C louds in the Aquarius simulations 
(|Springel et al.l[^QQ8f 1. Although their analysis ruled out 
Draco and UMi as possible members of the Magellanic 
group, they suggested that there could possibly still re¬ 
main a large number of faint satellites near the Clouds, 
in particular if the LMC and the SMC are on their first 
approach of the MW. 

There is also evidenc e of satellite associ ations around 
other nearby galaxies. Ilbata et ^ (j2Q13[ 1 argued that 
half of the satellites of M31 may lie in a thin, extended 
and rotating plane. Measurements of velocities of pairs 
of diametrically opposed satellite galaxies also hint that 
planes of satellites may be common in the low redshift 
Unive rse (jlbata et al.l Iml . Very recently, iTnllv et al.l 
(|2Q15f 1 has argued that all but 2 of the 29 galaxies in 
the Centaur us Group lie in one of two planes. Satel¬ 
lite galaxies are clearly not isotropically distributed and 
correlations are widespread. 

This suggests a picture in which some or all of the 
objects described in this paper - together with possibly 
some of the classical satellites like Draco and Ursa Minor 
- really did belong to a loose association with the LMC 
and the SMC as its largest members. The association 
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Figure 19. Distance to the Galactic center as a function of the 
distan ce t o the LMC. The symbols follow the same scheme as in 
Figure fTsI The blue-white 2D histogram in the background gives 
the expected density of objects assuming isothropic distribution on 
the sky (see text for more details). The red dashed line defines the 
’’zone of influence” of LMC. The inset shows the distribution of 
the number of objects inside the zone of influence as predicted by 
reshuffling the position vectors of the known satellites. The red 
solid line in the inset shows the actual number of satellites in the 
’’zone of influence” of LMC. We conclude that the detected objects 
constitute an overdensity around the LMC with a significance of 
94%. 


would then have a mass of 5 to 10 per cent of the mass of 
the Milky Way. Given the present-day distances, which 
range from 30 kpc (Reticulum 2) to 380 kpc (Eridanus 2) 
together with the bunching of seven satellites at 100 kpc, 
then the overall radius of the association must have been 
^150 kpc. These are very comparable to the masses and 
radii o f the dwarf associations discovered by iTullv et al.l 

(l2nnl . 

What might th e orbit of such a Magellanic Group be? 
iKallivavalil et al.l (j2013[ ) measured the proper motion of 
the LMC using a 7 year baseline. Orbit calculations de¬ 
pend on the uncertain masses of the Milky Way and the 
LMC, but the LMC is typically on first infall. Only if the 
mass of the Milky Way exceeds 2 x IO^^Mq and the LMC 
is comparatively light (3 x lOiOM©) can more than one 
pericentric passage occur. First infall solutions are also 
attractive as they also solv e the problem o f the anoma- 
lous blueness of the LMC. iTollerud et al.l (j2Qllf ) exam¬ 
ined a spectroscopic sample of isolated galaxies in SDSS, 
and found that bright satellite galaxies around Milky 
Way-type hosts are redder than field galaxies with the 
same luminsity. However, the LMC is unusually blue, 
which is plausibly a consequence of triggered star for¬ 
mation upon first infall. If the LMC is on first infall 
and if the mass of the Milky Way is 1 x IO^^Mq (e.g. 
iGibbons et al.l[2QT^ . then the orbital period is 9-12 Gyr. 
Increasing the mass of the Milky Way to 1.5 x lO^^M© 
gives orbital periods of 7-9 Gyr. This would imply that 
the Magellanic Group fell into the Milky Way halo at 
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Figure 20. Distribution of the MW satellites in the orbital plane 
of the LMC. The sym bols are assigned following the convention de¬ 
scribed in Figure flSl Coloured 2d histograms reflect the projection 
of the surveyed volume of SDSS (blue), VST ATLAS (green) and 
DES (lilac), with darker shades indicating nearly complete cover¬ 
age of the volume out to 250 kpc. Arrow shows the direc tion of 
the LMC’s motion as measured bv IKallivavalil et al.l (I2013I L Four 
LMC orbits are shown: two integrated forward in time and two 
backward in time for two different MW potential, one with a DM 
halo with concentration c = 10 and one with c = 24. 
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Figure 21. Distance from the LMC orbital plane \Zlmc\ as a 
function of the distance from the center of the LMC. Symbols are 
assigned as before, albeit the MW globular clusters are grey to aid 
clarity. Only three 3 of 9 new satellites show alignment with the 
LMC’s orbital plane. These are Reticulum 2, Horologium 1 and 
Eridanus 3. 

about a redshift oi z 1. 

Motivated by t he data of iTullv et akl (|2QQ6[) , 
iBowden et al ] (Iml looked at the survivability of loose 
associations in the outer parts of galaxies like the Milky 
Way. They found that the timescale required for an as¬ 
sociation to phase mix away is 10-15 times longer than 
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Figure 22. Positions o f the DES satel l ites w ith respect to the Magellanic gaseous Stream near the Southern Galactic Pole. The Stream’s 
HI column density from IPutman et aTl (|2003l) is shown as filled contours, with darker shades corresponding to higher densities. Hatched 
area contour shows the current DES footprint. The DES satellites appear to be avoiding regions with high HI column density. 


the orbital period (see their Figure 2). So, the struc¬ 
ture of such a Magellanic Group can therefore persist for 
enormous times, given the already long orbital periods. 

One prediction of the idea of a Magellanic Group is 
that the satellites might be expected to share the same 
sense of ciculation as the LMG and SMG about the 
Galaxy. However, on reflection, it is unclear whether this 
really need be the case. Modern simulations indicate that 
the LMG and SMC have recently come clos e enough to 
do spectacular daniage t o themselves (e.g., iBesla et ahl 
|2Q1Q[ iDiaz fc Bekkill2QlTI) The large-scale gaseous struc¬ 
tures of the Magellanic stream, the off-centered bar of 
the LMG and the irregular geometry of the SMC are 
all now believed to be the results of recent encounters 
between the LMG and the SMC, rather than any inter¬ 
action with the Milky Way. As such a violent interaction 
proceeds, any satellites or globular clusters of the LMG 
or SMC could be flung backwards against the prevailing 
circulation. 

A firmer - though rather less specific - prediction 
is that there must be more satellites in the vicinity of 
the Magellanic Clouds. We have shown that Reticulum 
2, Horologium 1, and Eridanus 3 are aligned with the 
LMC’s orbital plane and form part of the entourage of 
the LMG. Similarly, Tucana 2, Phoenix 2 and Grus 1 
most likely comprise part of the entourage of the SMC. 
However, all the objects in this paper trail the LMG and 
SMC. There must also be a counter-population of satel¬ 
lites that lead the LMG and SMC. Hunting is always eas¬ 
ier once we know where the big game is plentiful. And 
the happy hunting grounds for the ‘beasts of the south¬ 
ern wild’ are the Magellanic Clouds, especially in front 
of the LMC/SMC pair as reckoned by the LMG motion. 
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